On possible experimental realizations of directed percolation 
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Directed percolation is one of the most prominent universality classes of nonequilibrium phase tran- 
sitions and can be found in a large variety of models. Despite its theoretical success, no experiment 
is known which clearly reproduces the critical exponents of directed percolation. The present work 
compares suggested experiments and discusses possible reasons why the observation of the critical 
exponents of directed percolation is obscured or even impossible. 
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I. INTRODUCTION 

Physical phenomena far from thermal equilibrium are 
very common in nature. For example, many systems are 
subjected to an external flow of energy or particles which 
keeps them away from equilibrium. Similarly, dynamic 
systems starting with a nonequilibrium initial state may 
need a long time to reach thermal equilibrium. Theories 
of systems out of equilibrium are more difficult than equi- 
librium statistical mechanics since the partition sum is no 
longer given by the Gibbs ensemble. On the other hand, 
nonequilibrium systems may exhibit a potentially richer 



behavior than systems at thermal equilibrium. There- 
fore, the study of nonequilibrium phenomena is a field of 
growing interest, both theoretically and experimentally. 

A particularly interesting topic is the investigation of 
phase transitions far from equilibrium [n|. The theoret- 
ical interest in nonequilibrium phase transitions mainly 
stems from the emergence of universal features of the 
associated critical behavior. The concept of universal- 
ity was originally introduced by experimentalists in the 
context of equilibrium systems in order to describe the 
observation that order parameters of various apparently 
unrelated systems may display the same type of singu- 
lar behavior near the transition. These singularities are 
associated with a certain set of critical exponents which 
characterizes the universality class of the transition. The 
subsequent development of powerful theoretical concepts 
such as scaling, renormalization group, and conformal in- 
variance supported this hypothesis and established uni- 
versality as a paradigm of equilibrium statistical mechan- 
ics. 

Because of this success, theoretical physicists are nowa- 
days trying to transfer the idea of universality to nonequi- 
librium phase transitions. However, in the nonequilib- 
rium case the emerging picture remains less clear. Al- 
though universality certainly exists on the level of simple 
models, the experimental evidence of universal behav- 
ior under nonequilibrium conditions is still very poor. 
Therefore, it is not yet known to what extent the concept 
of universality can be applied to nonequilibrium critical 
phenomena. 

An important example is the universality class of Di- 
rected Percolation (DP) which describes continuous phase 
transitions from a spreading (wet) phase into an absorb- 
ing (dry) state B . The DP universality class is extremely 
robust with respect to the choice of the dynamic rules and 
covers a large variety of models with applications rang- 
ing from catalytic reactions g] and interface growth ^ to 
turbulence IpL The observed robustness led Janssen and 
Grassberger] g] to the conjecture that a continuous phase 
transition from a fluctuating active phase into a single ab- 
sorbing state should belong to the DP universality class, 
provided that the model uses short-range dynamics with- 
out special attributes such as additional symmetries or 



quenched randomness. In fact, DP is the canonical uni- 
versality class for nonequilibrium phase transitions into 
absorbing states. Thus it may be as important as the 
Ising universality class in equilibrium statistical mechan- 
ics. 

Despite this success in theoretical statistical physics, 
the critical behavior of DP, especially the set of critical 
exponents, has not yet been confirmed experimentally. 
The lack of experimental evidence is indeed surprising, 
especially since a large number of possible experimental 
realizations have been suggested in the past. As Grass- 
berger emphasizes in a summary on open problems in 
DP0: 

"...there is still no experiment where the crit- 
ical behavior of DP was seen. This is a very 
strange situation in view of the vast and suc- 
cessive theoretical efforts made to understand 
it. Designing and performing such an experi- 
ment has thus top priority in my list of open 
problems. ". 

The aim of the present work is to review the most im- 
portant experiments which have been suggested so far 
and to discuss their specific problems which could ob- 
scure the verification of the critical exponents. After 
a brief introduction to DP in Sect. || we will first dis- 
cuss certain catalytic reactions on two-dimensional sur- 
faces which mimic the dynamic rules of DP. As described 
in Sect. |^, DP may also be realized in certain wetting 
experiments where the interface between air and liquid 
undergoes a depinning transition. Another possible re- 
alization may be provided by systems of flowing sand 
on an inclined plane (see Sect. M). However, in none of 
these experiments the predicted critical exponents could 
be confirmed convincingly. 

What might be the reason for the apparent lack of ex- 
perimental evidence? It seems that the basic features of 
DP, which can easily be implemented on a computer, are 
quite difficult to realize in nature. One of these theoret- 
ical assumptions is the existence of an absorbing state. 
In real systems, however, a perfect non-fluctuating state 
cannot be realized. For example, a poisoned catalytic 
surface is not completely frozen, it will rather always be 
affected by small fluctuations. Although these fluctua- 
tions are strongly suppressed, they could still be strong 
enough to 'soften' the transition, making it impossible to 
quantify the critical exponents. 

Another reason might be the influence of quenched dis- 
order due to spatial or temporal inhomogeneities. In 
most experiments frozen randomness is expected to play 
a significant role. For example, a real catalytic surface is 
not fully homogeneous but characterized by certain de- 
fects leading to spatially quenched disorder. As shown in 
Ref. Q , this type of disorder may affect or even destroy 
the critical behavior of DP. Technical details concerning 
quenched disorder are summarized in the Appendix. 
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FIG. 1. Directed bond percolation on a diagonal square 
lattice. Open (closed) bonds are represented as solid (dashed) 
lines. The spreading agent, introduced at the site marked by 
the circle, flows downwards through open bonds, generating 
a certain cluster (bold bonds). 



II. DIRECTED PERCOLATION 

Directed percolation was introduced by Broadbent and 
Hammersley B in 1957 as an anisotropic variant of or- 
dinary percolation |10|] . The aim was to describe the 
connectivity of directed random media such as porous 
rocks in a gravitational field. Assuming that water can 
only move downwards, the permeability of a rock will de- 
pend on the average connectivity of its pores. One of the 
simplest models is directed bond percolation. As shown 
in Fig. m, neighboring sites (pores) of a diagonal square 
lattice are connected by bonds (channels) which are open 
with probability p and otherwise closed. Because of the 
gravitational field, the bonds function as 'valves' where- 
fore the spreading agent can only percolate along the 
given direction, as indicated by the arrows. 

Depending on the actual configuration of open bonds, 
each site generates a certain cluster of connected sites. 
A cluster of this kind would correspond to the maximal 
spreading range if the water was injected into a single 
pore. Below a certain threshold p < Pc all clusters are 
finite, i.e., the material is impermeable on large scales. 
However, above the critical value a cluster may become 
infinite so that water can percolate over arbitrarily long 
distances (cf. Fig||). 

Regarding the given direction as 'time', DP may be 
interpreted as a d+l-dimensional dynamic process de- 
scribing the spreading of some non-conserved agent. For 
example, we may enumerate the sites in Fig. m horizon- 
tally by a spatial coordinate i and vertically by a discrete 
time variable t. Then, for a given state at time t, we can 
determine the state at time i -I- 1 by means of certain 
stochastic updates pl| . 

Interpreting active sites as particles A and inactive 
sites as vacancies 0, the particles can either destroy 
themselves or produce an offspring. Moreover, if two 
particles reach the same site, they coagulate to a single 
particle. Therefore, DP may be regarded as a reaction- 
diffusion process 



diffusion: 

self-destruction: A 

offspring production: A 

coagulation: A 



A- 


^A 







A + A 


A- 


^A 







(1) 



position i 




P<Pc P=Pc P>Pc 

FIG. 2. Typical DP cluster starting from a single seed 
below, at, and above criticality. 

Depending on the ratio between offspring production and 
self destruction, the process may either remain active or 
reach the empty state from where it cannot escape. This 
is the so-called absorbing state of DP systems. For a spe- 
cial realization of DP, the so-called contact process [^ , 
the existence of a continuous transition between survival 
and extinction could be proven rigorously |59| ] . Near the 
transition the stationary density of active sites vanishes 
as a power law 



Pstat ~ (p - PcY 



(2) 



where /3 is the critical exponent associated with the order 
parameter p. Moreover, the DP process is characterized 
by a spatial and a temporal correlation length diverging 
at the transition as 



£.^^\p- Pc 



eii^b-pcl-"" 



(3) 



As already mentioned, a large variety of models show es- 
sentially the same properties at the transition, forming 
the DP universality class. The DP class corresponds to a 
specific field theory mM and is characterized by the three 
critical exponents /3, i/_l, i/||. Despite its simplicity, DP 
has not yet been solved exactly. However, the critical 
exponents can easily be estimated by computer simula- 
tions. The most accurate estimates are summarized in 
Table L Notice that in d > 4 spatial dimensions fluctua- 
tions become irrelevant so that the exponents are given 
by their mean field values. 



III. CATALYTIC REACTIONS 

It is well known that under specific conditions cer- 
tain catalytic reactions mimic the microscopic rules of 
DP models. For example, the Ziff-Gulari-Barshad (ZGB) 
model, which was designed in order to describe the cat- 
alytic reaction CO -I- O ^ CO2 on a platinum surface Q, 
displays a DP transition. In the ZGB model a gas com- 
posed of CO and O2 molecules with fixed concentrations 
y and I — y, respectively, is brought into contact with a 
catalytic material. The catalytic surface is represented 
by a square lattice whose sites can either be vacant (0), 
occupied by a CO molecule, or occupied by an O atom. 
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TABLE I. Numerical estimates for the critical exponents 
of directed percolation in d+1 dimensions. 



CO molecules fill any vacant site at rate y, whereas O2 
molecules dissociate on the surface into two O atoms and 
fill pairs of adjacent vacant sites at rate I — y. Finally, 
neighboring CO molecules and O atoms recombine in- 
stantaneously to CO2 and desorb from the surface. On 
the lattice the three processes correspond to the reaction 
scheme 
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(4) 



Therefore, if the whole lattice is entirely covered either 
with CO or O, the system is trapped in a catalytically in- 
active state. These 'poisoned' states are the two absorb- 
ing configurations of the ZGB model. As shown in Fig. ||, 
the corresponding phase diagram displays two absorbing 
phases. For y < yi ~ 0.389 the system evolves into 
the 0-poisoned state whereas for y > y2 — 0.525 it al- 
ways reaches the C 0-poisoned state. Between these two 
values the model is catalytically active. The two transi- 
tions into the absorbing phases are different in character, 
namely discontinuous at y = y2 and continuous aX y = yi 
(see Fig, m. Motivated by the DP conjecture, Grinstein 
et al. ||l^ expected the latter to belong to the DP univer- 
sality class. In order to verify this hypothesis, extensive 
numerical simulations were performed. Initially it was 
believed that the critical exponents were different from 
those of DP ||l^ , while later the transition at j/ = yi was 
found to belong to DP ||l^. Very precise estimates of 
the critical exponents were recently obtained in Ref . p8| , 
confirming the DP conjecture. DP exponents were also 
obtained in a simplified version of the ZGB model p9[ . 
These theoretical models therefore suggest that certain 
catalytic reactions could serve as an experimental real- 
ization of DP. 

In real catalytic reactions, however, only the discontin- 
uous transition aXy — y2 can be observed. The schematic 
graph on the right-hand side of Fig. g shows the reaction 
rates as functions of the CO pressure measured in a cat- 
alytic reaction on a Pt(210) surface [g2[. Although the 
experiment was designed in order to investigate the tech- 
nologically interesting regime of high activity close to the 
first-order phase transition, it clearly indicates that poi- 
soning with oxygen does not occur. Instead the reactivity 
increases almost linearly with the CO pressure. Similar 
results were obtained for Pt(lll) and for other catalytic 
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FIG. 3. Catalytic reactions in theory (left) and experi- 
ment (right). The schematic graphs show the concentrations 
of oxygen (solid line) and carbonmonoxide (dashed line) as a 
function of the CO adsorption rate. 



materials. Thus, so far there is no experimental evidence 
for the DP transition predicted by the ZGB model. 

One may speculate why the DP transition is obscured 
or even destroyed under experimental conditions. A pos- 
sible reason might be the reaction chain being much more 
complicated than in the ZGB model |2^]. Moreover, the 
0-poisoned system might not be in a perfect absorb- 
ing state, i.e., the surface can still adsorb CO molecules 
although it is satured with carbonmonoxide. Another 
possibility is thermal (nonreactive) desorption of oxy- 
gen which - in the DP language - would correspond to 
spontaneous creation of active sites due to an external 
field p3]. Finally, defects and inhomogeneities of the 
catalytic material could lead to an eiTective (spatially 
quenched) disorder. As shown in the Appendix, this type 
of disorder is marginal. It can therefore seriously affect 
the critical behavior and even modify the values of the 
exponents. 

For a long time microscopic details were difficult to 
study experimentally. However, novel techniques such as 
scanning tunneling microscopy (STM) led to an enor- 
mous progress in the understanding of catalytic reac- 
tions. They also point at various unexpected subtleties. 
For example, recent experiments revealed that the reac- 
tions preferably take place at the perimeter of oxygen 
islands ||2^. Furthermore, it was observed that the ad- 
sorbed CO molecules on Pt(lll) may form three different 
rotational patterns representing the c(4x2) structure of 
CO on platinum, leading to three competing absorbing 
states |26|. Moreover, the STM technique allows one to 
trace individual molecular reactions and to determine the 
corresponding reaction rates. In addition, the influence 
of defects such as terraces on catalytic reactions can be 
quantified experimentally p7[ . We may therefore expect 
a considerable progress in the understanding of catalytic 
reactions in near future. 



IV. GROWING INTERFACES 

Various models for interface growth exhibit a roughen- 
ing or depinning transition which can be related to DP. 
In this Section we discuss three examples, namely tran- 
sitions of depinning interfaces in random media, polynu- 
clear growth processes, and solid-on-solid growth with 
evaporation at the edges of plateaus. 



A. Depinning transitions 

Depinning transitions of driven interfaces provide a 
very promising class of experiments which could be re- 
lated to DP ||2^ . In these experiments a liquid is pumped 
through a porous medium. If the driving force F is suffi- 
ciently low the liquid cannot move through the medium 
since the air/liquid interface is pinned at certain pores. 
Above a critical threshold, however, the interface starts 
moving through the medium with an average velocity v. 
Close to the transition, v is expected to scale as 



F-F, 



(5) 



where 9 is the velocity exponent. Moreover, in the mov- 
ing phase F > Fc the interface roughness averaged over 
length £ should obey the usual scaling law for roughening 
interfaces |E9|l 



w{i,t)^rf{t/n, 



(6) 



where a is the so-called roughening exponent. One of 
the first experiments in l-l-I dimensions was performed 
by Buldyrev et al., who studied the wetting of paper 
in a basin filled with suspensions of ink or coffee pOJ. 
Measuring the interface width they found the roughness 
exponent a — 0.63(4). In various other experiments the 
values are scattered between 0.6 and 1.25. This is surpris- 
ing since the Kardar-Parisi-Zhang (KPZ) class pl| , p2[ , 
the canonical universality class for roughening interfaces, 
predicts the exponent a = 1/2 which is smaller than the 
experimentally observed values. 

It is believed that the large values of a are due to in- 
homogeneities of the porous medium. Due to these inho- 
mogeneities, the interface does not propagate uniformly 
by local fluctuations as in the KPZ equation, it rather 
propagates by avalanches. In the literature two univer- 
sality classes for this type of interfacial growth have been 
proposed. In case of linear growth the interface should 
be described by a random field Ising model ||3^ , leading 
to the exponents a = 1, z = 4/3, and 9 — 1/3 in 1-1-1 
dimensions. In the presence of a KPZ-type nonlinearity, 
however, the roughening process should exhibit a depin- 
ning transition which is related to DP [[4|j34[| . Notice that 
the underlying DP mechanism of depinning transitions 




FIG. 4. Simple model exhibiting a depinning transition. 
The pores are represented by cells on a diagonal square lattice. 
The permeability across the edges of the cells depends on the 
direction of flow: In downwards direction all edges are per- 
meable whereas in upwards direction they are permeable with 
probability p and impermeable otherwise. The right panel of 
the figure shows a particular configuration of open (dashed) 
and closed (solid) edges. Pumping in water from below, the 
interface becomes pinned along a eirected path of solid lines, 
leading to a finite cluster of wet cells (shaded region). The 
dashed arrow represents an open path in order to illustrate 
the flow. 



differs significantly from an ordinary directed percolation 
process in a porous medium subjected to a gravitational 
field (cf. Sect. VI). In the latter case the spreading agent 
is restricted to percolate along a given direction, i.e., the 
flow is strictly unidirectional. In depinning experiments, 
however, water may flow both according the pumping 
force and - even more easily - in opposite direction, as 
will be explained below. 

A simple model exhibiting a depinning transition is 
shown in Fig. ^. In this model the pores are represented 
by cells of a diagonal square lattice. The liquid can flow 
to neighboring cells by crossing the edges of the cell. De- 
pending on the direction of the flow these edges can ei- 
ther be open or closed. For simplicity we assume that all 
edges are permeable in downwards direction, whereas in 
upwards direction they can only be crossed with a certain 
probability p. Thus, by starting with a horizontal row of 
wet cells at the bottom, we obtain a compact cluster of 
wet cells, as illustrated in Fig. H The size of the cluster 
(and therewith the penetration depth of the liquid) de- 
pends on p. If p is large enough, the cluster is infinite, 
corresponding to a moving interface. If p is sufficiently 
small, the cluster is bound from above, i.e., the interface 
becomes pinned. 

The depinning transition is related to DP as follows. 
As can seen in the figure, a pinned interface may be 
interpreted as a directed path along impermeable edges 
running from one boundary of the system to the other. 
Obviously, the interface becomes pinned only if there ex- 
ists a directed path of impermeable bonds connecting the 
boundaries of the system. Hence the depinning transition 
is related to an underlying DP process running perpen- 
dicular to the direction of growth. The pinning mech- 
anism is illustrated in Fig. |g, where a supercritical DP 
cluster propagates from left to right. The cluster's back- 
bone, consisting bonds connecting the two boundaries, is 




FIG. ■"). Pimietl iutcifacc and the uudorlyiiig DP pKjccsf 
The figure is explained in the text. 



indicated by bold dots. The shaded region denotes the 
resulting cluster of wet cells. As can be seen, the inter- 
face will be pinned at the lowest lying branch of the DP 
backbone. Therefore, the roughening exponent coincides 
with the meandering exponent 



t^J-/i^|| 



(7) 



of the backbone. Moreover, by analyzing the dynanfics of 
the moving interface, it can be shown that the dynamic 
critical exponents are given hy 9 — a and z — 1. Thus, 
depinning transitions in inhomogeneous porous media 
may serve as experimental realizations of the DP uni- 
versality class. 

Comparing the prediction (|7|) with the result a — 
0.63(4) obtained by Buldyrev et al. Q we find an ex- 
cellent coincidence, confirming the validity of the model 
introduced above. Therefore, this experiment can be re- 
garded as a first experimental evidence of DP exponents. 
However, only one exponent has been verified, and it is 
not fully clear how accurate and reproducible these ex- 
ponents are. Further experimental effort in this direction 
would be desirable. 

Similar experiments were carried out in 2+1 dimen- 
sions with a spongy-like material used by florists, as well 
as fine-grained paper rolls |^. In this case, however, the 
exponent a is not related to 2+1-dimensional DP, instead 
it corresponds to the dynamic exponent of percolating di- 
rected interfaces in 2+1 dimensions. In experiments as 
well as in numerical simulations a roughness exponent 
a — 0.50(5) was obtained. 

In order to perform a depinning experiment which can 
easily be reproduced, Dougherty and Carle measured the 
dynamical avalanche distribution of an air/water inter- 
face moving through a porous medium made of glass 
beads |3g]. Assuming an underlying DP process, the dis- 
tribution P{s) of avalanche sizes s is predicted to behave 
algebraically. In the experiment, however, a stretched 
exponential behavior P{s) ~ s~^e~*/^ is observed even 
for small fiowing rates. The estimates for the exponent 
b are inconclusive; they depend on the time window of 
the measurement and vary between —0.5 and 0.85. Even 
more recently Albert et al. proposed a method allowing 
to identify the universality class by measuring the propa- 
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FIG. 6. Polynuclear growth model. In the first half time 
step atoms are deposited with probability p. In the second 
half time step islands grow deterministically by one step and 
coalesce. 



gation velocity of locally tilted parts of the interface [^ . 
Their results suggest that interfaces propagating in glass 
beads are not described by a DP depinning process, but 
to be related instead to the random-field Ising model. 
From the theoretical point of view this is surprising since 
linear growth of the interface is a special case which re- 
quires fine-tuning of certain parameters. Further experi- 
mental effort would be needed to understand these find- 
ings. 



B. Polynuclear growth 

A completely different DP mechanism is responsible for 
roughening transitions in so-called polynuclear growth 
(PNG) models [||-|o|. A key feature of PNG models 
is the use of parallel updates, leading to a maximal prop- 
agation velocity of one monolayer per time step. For a 
high adsorption rate the interface is smooth, propagat- 
ing at maximal velocity v — I. Decreasing the adsorp- 
tion rate below a certain critical threshold, PNG models 
exhibit a roughening transition to a rough phase with 
V < 1. In contrast to equilibrium roughening transitions, 
which only exist in d > 2 dimensions, PNG models have 
a roughening transition even in one spatial dimension. 

Probably the simplest PNG model investigated so far 
is defined by the following dynamic rules [p8| . In the 
first half time step atoms 'nucleate' stochastically at the 
surface by 



h,{t+ 1/2) = 



hi{t) + 1 with prob. p , 
hi{t) with prob. 1— p 



(8) 



In the second half time step the islands grow determin- 
istically in lateral direction by one step. This type of 
growth may be expressed by the update rule 

h,{t + l)= max [/ij(i-M/2),/ij(t-M/2)l , (9) 

where j runs over the nearest neighbors of site i. 

The relation to DP can be established as follows. Start- 
ing from a flat interface ft.,;(0) = 0, let us interpret sites 
at maximal height hi{t) — t as active sites of a DP pro- 
cess. The adsorption process (H) turns active sites into 
the inactive state with probability 1 — p, while the pro- 
cess (0) resembles offspring production. Therefore, if p 
is large enough, the interface is smooth and propagates 



with maximal velocity v = 1. This situation corresponds 
to the active phase of DP. Therefore, we expect the den- 
sity of sites at maximal height to scale as 



^i:*' 



-t^ [p- PcY 



(10) 



where N denotes the system size. Below a critical thresh- 
old, however, the density of active sites at the maximum 
height hi{t) = t vanishes; the growth velocity is smaller 
than 1 and the interface evolves into a rough state. Al- 
though this mapping to DP is not exact, numerical simu- 
lations suggest that it still remains valid on a qualitative 
level. More specifically, it turns out that PNG models 
are related to a unidirectionally coupled hierarchy of DP 
processes |pl| . 

A closely related class of models was introduced in or- 
der to describe the growth of colonial organisms such 
as fungi and bacteria [|42| , motivated by recent experi- 
ments with the yeast Pichia membranaefaciens on solid- 
ified agarose film |^^. By varying the concentration of 
polluting metabolites, different front morphologies were 
observed. The model proposed in |^^ aims to explain 
these morphological transitions on a qualitative level. A 
careful analysis of the dynamic rules shows that mod- 
els for fungal growth and PNG models are very similar 
in character. They both employ parallel dynamics and 
exhibit a DP-related roughening transition. 

Concerning experimental realizations of PNG models, 
one major problem ~ apart from quenched disorder - is 
the use of parallel updates. The type of update is cru- 
cial; by using random-sequential updates the transition 
is lost since in this case there is no maximum velocity. 
However, in realistic experiments atoms do not move 
synchronously, but the adsorption events are rather ran- 
domly distributed in time. Therefore, random sequential 
updates might be more appropriate to describe such ex- 
periments. It thus remains an open question to what 
extent PNG processes can be realized in nature. In fact, 
it would be interesting to see if PNG can be generalized 
to random-sequential dynamics. 



C. Growth with evaporation at the edges of plateaus 

DP-related roughening transitions can also be observed 
in certain solid-on-sol id g rowth processes with random- 
sequential updates Q,^ . As a key feature of these mod- 
els, atoms may desorb exclusively at the edges of exist- 
ing layers, i.e., at sites which have at least one neigh- 
bor at a lower height. By varying the growth rate, such 
growth processes display a roughening transition from a 
non-moving smooth phase to a moving rough phase. 

A simple solid-on-solid model for this type of growth 
is defined by the following dynamic rules Q]: For each 
update a site i is chosen at random and an atom is ad- 
sorbed 
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FIG. 7. Restricted solid-on-solid growth model exhibiting 
a roughening transition from a non-moving smooth to a mov- 
ing rough phase. Monomers are randomly deposited whereas 
desorption takes place only at the edges of plateaus. 



hi —> hi + 1 with probability q 

or desorbed at the edge of a plateau 

hi —f imn{hi, /li+i) with probability (1 — q)/2 
hi ~f imn{hi, /ii_i) with probability (1 — q)/2 



(11) 



(12) 



Moreover, the growth process is assumed to be restricted, 
i.e., updates are only carried out if the resulting configu- 
ration obeys the constraint 



ti±i 



< 1 



(13) 



The qualitative behavior of this model is illustrated in 
Fig. 0. For small q the desorption processes ( |l2| ) dom- 
inate. If all heights are initially set to the same value, 
this level will remain the bottom layer of the interface. 
Small islands will grow on top of the bottom layer but will 
be quickly eliminated by desorption at the island edges. 
Thus, the interface is effectively anchored to its bottom 
layer and a smooth phase is maintained. The growth ve- 
locity V is therefore zero in the thermodynamic limit. As 
q is increased, more islands on top of the bottom layer 
are produced until above qc ~ 0.189, the critical value 
of q, they merge forming new layers at a finite rate, giv- 
ing rise to a finite growth velocity. In an infinite system 
the growth velocity scales near the transition as 



^i'-{p~Pcr^ 



(14) 



Since the propagation velocity fluctuates locally, the in- 
terface evolves into a rough state according to the pre- 
dictions of the KPZ equation. The same type of critical 
behavior is observed in similar models without the re- 
striction (p^). 

At the transition the dynamics of the model is related 
to DP as follows. Starting with a flat interface at zero 
height, let us consider all sites with hi = as parti- 
cles A of a DP process. Growth according to Eq. ( |ll| ) 
corresponds to spontaneous annihilation A — > 0. Con- 
versely, desorption may be considered as a particle cre- 
ation process. However, since atoms may only desorb 
at the edges of plateaus, particle creation requires a 




FIG. 8. Possible experimental realization of the growth 
process defined in Eqs. (fni) and (113) by lateral growth of a 
monolayer on a substrate. 



neighboring active site, corresponding to offspring pro- 
duction A -^ 2A. These rules resemble (although not 
exactly) the dynamics of a DP process. In contrast to 
PNG models, the DP process takes place at the bottom 
layer of the interface. Moreover, the roughening tran- 
sition does not depend on the use of either parallel or 
random-sequential updates. However, if composite par- 
ticles instead of monomers are deposited on the surface, 
the universality class may change due to additional sym- 
metries [{l6| . 

With respect to experimental realizations of the dy- 
namic rules ([ll])-(p^) we note that atoms are not allowed 
to diffuse on the surface. This assumption is rather un- 
natural since in most experiments the rate for surface 
diffusion is much higher than the rate for desorption 
back into the gas phase. Therefore, it will be difflcult 
to realize this type of homoepitaxial growth experimen- 
tally. However, in a different setup, the above model 
could well be relevant |4^. As illustrated in Fig. g, a 
laterally growing monolayer could resemble the dynamic 
rules (|ll| ) and (12) by identifying the edge of the mono- 
layer with the interface of the growth model. In this 
case 'surface diffusion', i.e. diffusion of atoms along the 
edge of the monolayer, is highly suppressed. Moreover, 
in single-step systems (such as fcc(lOO) surfaces) it would 
also be possible to implement the restriction (p^. 



V. FLOWING GRANULAR MATTER 

It has been shown recently that simple systems of flow- 
ing sand on an inclined plane, e^. the experiments per- 
formed by Douady and Daerr P8|,y , could serve as ex- 
perimental realizations of DP pO| . In the Douady-Daerr 
experiment glass beads with a diameter of 250-425 /im 
are poured uniformly at the top of an inclined plane cov- 
ered by a rough velvet cloth (see Fig. g). As the beads 
flow down, a thin layer settles and remains immobile. 
Increasing the angle of inclination by Acj) the layer be- 
comes dynamically unstable, i.e., by locally perturbing 
the system at the top of the plane an avalanche of flow- 
ing granular matter will be generated. In the experiment 
these avalanches have the shape of a fairly regular tri- 
angle with an opening angle 9. As the increment Aip 
decreases, the value of 6 decreases, vanishing as 




CDP: 



DP: 



FIG. 9. Simplified drawing of tlie Douady-Daerr experi- 
ment. At a given angle (j> a layer of a certain thickness settles. 
Perturbing the layer locally with a stick leads to an avalanche 
of flowing sand. 



tan 61 - (A(p)^ 



(15) 



with a certain critical exponent x. The experimental re- 
sults suggest the value x —1 p^ . 

In order to explain the experimentally observed trian- 
gular form of the avalanches, Bouchaud et al. proposed a 
mean-field theory based on deterministic equations tak- 
ing the actual local thickness of the flowing avalanche 
into account JSIJ. This theory predicts the exponent 
X = 1/2. Another explanation assumes that flowing sand 
may be associated with a nearest-neighbor spreading pro- 
cess pG] . Considering the avalanche as a cluster of active 
sites while identifying the vertical coordinate of the plane 
with time and the increment of inclination /S.if with p—pc 
the opening angle is expected to scale as 



tan 



e ^ e±/Cii ~ (A^)''" 



(16) 



where vu and iy± are the scaling exponents of the spread- 
ing process under consideration. 

To support this scaling argument, a simple lattice 
model was introduced which mimics the physics of flow- 
ing sand ||50|| . The model exhibits a transition from an in- 
active to an active phase with avalanches whose compact 
shapes reproduce the experimental observations. On lab- 
oratory scales the model predicts the exponent a; = 1, 
corresponding to the universality class of compact di- 
rected percolation (CDP) which is characterized by the 
exponents p2[ 



/3 = 0. 



(17) 



The CDP behavior, however, is only an initial transient 
and crosses over to DP after a very long time. Thus 
the Douady-Daerr experiment - performed on sufficiently 
large scales - may serve as a physical realization of DP. Ir- 
regularities of the layers thickness may affect the spread- 
ing properties of avalanches. However, such inhomo- 
geneities can be considered as a kind of spatio- temporally 
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FIG. 10. Typical clusters generated at criticality on small 
and large scales, illustrating the crossover from CDP to DP. 

quenched disorder which is irrelevant on large scales (see 
Appendix). Thus, in contrast to the previous examples, 
the problem of quenched disorder does not play a major 
role in this type of experiments. 

The crossover from CDP to DP is very slow and 
presently not accessible in the experiments. To illus- 
trate the crossover, two avalanches are plotted on dif- 
ferent scales in Fig. nQ. The left one represents a typical 
avalanche within the first few thousand time steps. As 
can be seen, the cluster appears to be compact. However, 
as shown in the right panel of Fig. |l^, the cluster breaks 
up into several branches after a very long time. As a pre- 
condition for DP behavior, initially compact avalanches 
should thus be able to break up into several branches. 
Only then is it worthwhile to optimize the experimental 
setup and to measure the critical exponents quantita- 
tively. 

More recent experimental studies ||5^ confirm that 
for high angles of inclination critical avalanches do split 
up into several branches (see Fig. ^l]). Yet here the 
avalanches have no well defined front, the propagation ve- 
locity of separate branches rather depends on their thick- 
ness. It is therefore no longer possible to interpret the 
vertical axis as a time coordinate. Moreover, it is not 
yet known how the spreading process depends on corre- 
lations in the initial state. As shown in Ref. pM, such 
long-range correlations may change the values of certain 
dynamic critical exponents. However, recent studies of a 
single rolling grain on an inclined rough plane |55|| sup- 
port that there are presumably no long-range correlations 
due to a 'memory' of rolling grains. By means of molec- 
ular dynamics simulations it was shown that the motion 
of a rolling grain consists of many small bounces on each 
grain of the supporting layer. Therefore, the rolling grain 
quickly dissipates almost all of the energy gain from the 
previous step and thus forgets its history very fast. For 
this reason it seems to be unlikely that quenched disor- 
der of the prepared layer involves long-range correlations. 
Therefore, fiowing granular matter seems to be a promis- 
ing candidate for an experimental realization of DP. 
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FIG. 11. Avalanches splitting up into several branches, 
observed in recent experiments with high angle ol inclination 
(reprinted with kind permission Irom A. Daerr). 



VI. OTHER APPLICATIONS 

This section discusses other applications which are less 
promising to serve as experimental realizations of DP, 
although they are frequently quoted in the literature. 

Porous media: One example is percolating water in 
a porous medium subjected to an external driving force. 
The medium could be a porous rock in a gravitational 
field where neighboring pores are connected by channels 
with a certain probability. Depending on this probabil- 
ity, the penetration depth is either finite or the water 
may "percolate" over infinitely long distances through 
the medium. Due to the external driving force, the flow 
in the medium is assumed to be strictly unidirectional, 
i.e., the water can only flow downwards (in contrast to 
the depinning models of Sect. EVA). Although this appli- 
cation is quite natural, it is extremely difficult to realize 
experimentally. For example, by studying natural sand- 
stone a broad distribution of pore sizes was observed [ p6[ . 
Although these experiments are concerned with isotropic 
percolation, similar difficulties are expected in the di- 
rected case. 



Epidemics: Another frequently quoted applications 
of DP is the spreading of epidemics without immuniza- 
tion ||57|-|59|]. Here infection and recovery resemble the 
reaction-diffusion scheme (p. If the rate of infection is 
very low, the infectious disease will disappear after some 
time. If infections occur more frequently, the disease 
may spread and survive for a very long time. However, 
spreading processes in nature are usually not homoge- 
neous enough to reproduce the critical behavior of DP. 
Moreover, in many realistic spreading processes short- 
range interactions are no longer appropriate. This situa- 
tion emerges, for example, when an infectious disease is 
transported by insects. The motion of the insects is typ- 
ically not a random walk, rather occasional fiights over 
long distances may occur before the next infection takes 
place. On the level of theoretical models such long-range 
interactions may be described by Levy flights |p3] , lead- 



ing to continuously varying critical exponents | |63| . 

Forest fires: A closely related problem is the spread- 
ing of forest fires [ pO] . Tephany et al. studied the propa- 
gation of fiame fronts on a random lattice both under qui- 
escent conditions and in a wind tunnel pU . The experi- 
mental estimates of the critical exponents at the spread- 
ing transition are in rough agreement with the predic- 
tions of isotropic and directed percolation, respectively. 
However, the accuracy of these costly experiments re- 
mains limited. 

Calcium dynamics: DP transitions may also occur 
in certain kinetic models for the dynamics of Calcium 
ions in living cells. Ca^^ ions play an important phys- 
iological role as second messenger for various purposes 
ranging from hormonal release to the activation of egg 
cells by fertilization 1 64 6^ . The cell uses nonlinear prop- 



agation of increasing intracellular Ca^+ concentration, 
a so-called calcium wave, as a tool to transmit signals 
over distances which are much longer than the diffusion 
length. For example, propagating Ca^+ waves can be 
observed in the immature Xenopus laevis oocyte [p6|. 
So far theoretical work focused mainly on determinis- 
tic reaction-diffusion equations in the continuum which 
explain various phenomena such as solitary and spiral 
waves p7[ |. This mean-field type approach, however, ig- 
nores the infiuence of fluctuations. Yet, near the tran- 
sition between survival and extinction of Ca^"*", activity 
fluctuations may play an important role. Recently im- 
proved models have been introduced which take also the 
stochastic nature of Calcium release into account |£8|,^ . 
As expected, the transition in one of these models be- 
longs to the DP universality class p3]. However, from 
the experimental point of view it seems to be impossi- 
ble to confirm or disprove this conjecture. On the one 
hand, the size of a living cell is only a few order of mag- 
nitude larger than the diffusion length, leading to strong 
flnite-size effects in the experiment. On the other hand, 
inhomogeneities as well as internal structures of the cell 
give rise to a completely unpredictable form of quenched 
noise which may be correlated in space and time. There- 
fore, it seems to be impossible to identify the universality 
class of the transition in actual experiments. It would be 
rather an achievement to find clear evidence for the very 
existence of a phase transition between survival and ex- 
tinction of propagating calcium waves. For quantitative 
experiments, it would be interesting to reproduce the dy- 
namics of Calcium in a well-defined environment, e.g., on 
an artificial membrane |24[ . 

Directed polymers: DP is also related to the prob- 
lem of directed polymers U^ . In contrast to DP, which 
is defined as a local process, the directed polymer prob- 
lem selects directed paths in a random medium by global 
optimization. Under certain conditions, namely a bi- 
modal distribution of random numbers, both problems 
were shown to be closely related [Q. More specifically, 
the roughness exponent of the optimal path in a directed 
polymer problem is predicted to cross over from the KPZ 



value 2/3 to the DP value v^jv^ ~ 0.63 at the transi- 
tion point. Directed polymers were used to describe the 
propagation of cracks ¥%■ However, it is rather unlikely 
that crack experiments can reproduce the tiny crossover 
from KPZ to DP. 

Turbulence: Finally, DP has also been considered as 
a toy model for turbulence. According to Ref. B, the 
front between turbulent and laminar flow should exhibit 
the critical behavior of DP. For example, the velocity of 
the front should scale algebraically with a combination 
of DP exponents. However, these predictions are based 
rather on heuristic arguments than on rigorous results. 
In fact, in many respects turbulent phenomena show a 
much richer behavior than DP. 



VII. CONCLUSIONS 

Directed percolation has kept theoretical physicists fas- 
cinated since more than forty years. Several reasons make 
directed percolation so appealing. First of all, DP is a 
very simple model in terms of its dynamic rules. Never- 
theless, the DP phase transition turns out to be highly 
nontrivial. In fact, DP belongs to the very few critical 
phenomena which have not yet been solved exactly in one 
spatial dimension. Therefore, the critical exponents are 
not yet known analytically. High-precision estimates in- 
dicate that they might be given rather by irrational than 
by simple fractional values. 

Moreover, DP is extremely robust. It stands for a 
whole universality class of phase transitions from a fluc- 
tuating phase into an absorbing state. In fact, a large 
variety of models displays phase transitions belonging to 
the DP universality class. Thus, on the theoretical level, 
DP plays the role of a standard universality class similar 
to the Ising model in equilibrium statistical mechanics. 

In spite of its simplicity, no experiment is known con- 
firming the values of the critical exponents quantitatively. 
An exception may be the wetting experiment performed 
by Buldrey et al. where the value of the roughness expo- 
nent OL coincides with v^jv^ within less than 10%. How- 
ever, since the results of similar experiments are scattered 
over a wide range, further experimental effort would be 
needed in order to confirm the existence of DP in this 
type of systems. 

Apart from difficulties to realize a non-fluctuating ab- 
sorbing state, a fundamental problem of DP experiments 
is the emergence of quenched disorder due to certain in- 
homogeneities of the system. Depending on the type 
of disorder, even weak inhomogeneities might obscure 
or even destroy the DP transition. Therefore, the most 
promising experiments are those where quenched disor- 
der is irrelevant on large scales. This is th e case, for ex- 
ample, in wetting experiments (Sect. |lV A ) and systems 
of flowing granular matter (Sect. M). 

In spite of all these difficulties, many physicists believe 



that DP should have a counterpart in reality, mostly be- 
cause of its simplicity and robustness. Therefore, Grass- 
bergers message remains valid: The experimental real- 
ization of DP is an outstanding problem of top priority. 
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APPENDIX A: QUENCHED DISORDER 

On a coarse-grained scale the temporal evolution of a 
DP process without quenched disorder is described by 
the Langevin equation [|6| 



dtp = ap- Xp^ 



DV'^p 



VpC, 



(Al) 



where ^(x, t) denotes an uncorrelated Gaussian noise: 

(C(x, t)C(x', t')) = r<5'^(x - x')<5(i - t') . (A2) 

The noise C(x,t) represents the intrinsic fluctuations of 
the DP process due to the stochastic nature of the dy- 
namic rules. The parameter a controls the rate for off- 
spring production and can be thought of as a being mea- 
sure of the percolation probability p — Pc- 

Quenched disorder may be introduced by random vari- 
ations of the parameter a, i.e., by adding another noise 
field T]: 



a + T] . 



(A3) 



Thus, the resulting Langevin equation reads 

dtp = ap- Ap2 + dV^p -f VpC + PV ■ (A4) 

The noise rj is quenched in the sense that quantities like 
the particle density are averaged over many independent 
realizations of the intrinsic noise C, while the disorder field 
rj is kept fixed. In the following we distinguish three 
different types of quenched disorder: 

A. Spatially quenched disorder »ys(x). 

B. Temporally quenched disorder rjtit). 

C. Spatio-temporally quenched disorder /^^^(x, i). 

These variants of quenched disorder differ in how far they 
affect the critical behavior of DP. In the following we 
review some of the main results. 
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A. Spatially quenched disorder 

For spatially quenched disorder, the noise field rj is 
defined through the correlations 



r;,(x)7?,(x') = 7'5''(x-x'), 



(A5) 



where the bar denotes the average over independent re- 
alizations of the disorder field (in contrast to averages 
(. . .) over the intrinsic noise (,). The parameter 7 is an 
amplitude which controls the intensity of disorder. In 
order to find out whether this type of noise affects the 
critical behavior of DP, let us consider the properties of 
the Langevin equation under the scaling transformation 



Aj 



t^AH, 



A- 



V, 



(A6) 



where A is a dilatation factor while z = i^\\/i'± and 
X — P/v± are the critical exponents of DP. In absence of 
quenched disorder, the Langevin equation turns out to 
be invariant under rescaling if 



^ = 2, 



X = 2, 



a = 0. 



d — dr 



(A7) 



where dc = 4 is the upper critical dimension of DP. These 
values are consistent with the DP mean field exponents 
(3 — \, v^_ ~ 1/2, and v^ = 1, which are valid in d > 4 
dimensions. Checking the scaling behavior of the addi- 
tional term pr]g in Eq. (A4) at the critical dimension, we 
observe that it scales as 



PVs 



\-da/2-x 



PVs, 



(A8) 



i.e., spatially quenched disorder is a marginal perturba- 
tion. Therefore, it may seriously affect the critical be- 
havior at the transition. 

The same result is obtained by considering the field- 
theoretic action. Without quenched noise, DP is de- 
scribed by the action of Reggeon field theory pj] 



So^ I d'^x I dti; dt 



DV^ + g{iP - iP) V (A9) 



where '0(x, i) represents the local particle density while 
'0(x, i) denotes the Martin-Siggia- Rosen response field. 
As shown by Janssen |7^], spatially quenched noise can 
be taken into account by adding the term 



S = So + j 



d'^x 



dtipip 



(AlO) 



By simple power counting we can prove that this addi- 
tional term is indeed a marginal perturbation. Janssen 
showed by a field-theoretic analysis that the stable fixed 
point is shifted to an unphysical region, leading to run- 
away solutions of the flow equations in the physical re- 
gion of interest. Therefore, spatially quenched disor- 
der is expected to crucially disturb the critical behav- 
ior of DP. The findings are in agreement with earlier 



numerical results by Moreira and Dickman [Q who 
reported non-universal logarithmic behavior instead of 
power laws. Later Cafiero et al. [Q showed that DP with 
spatially quenched randomness can be mapped onto a 
non-Markovian spreading process with memory, in agree- 
ment with previous results. 

From a more physical point of view, spatially quenched 
disorder in 1-1-1 dimensional systems was studied by Web- 
man et al. [g[ . It turns out that even very weak random- 
ness drastically modifies the phase diagram. Instead of 
a single critical point one obtains a whole phase of very 
slow glassy-like dynamics. The glassy phase is charac- 
terized by non-universal exponents which depend on the 
percolation probability and the disorder amplitude. For 
example, in a supercritical 1-1-1 dimensional DP process 
without quenched disorder the boundaries of a cluster 
propagate at constant velocity v. However, in the glassy 
phase V decays algebraically with time. The correspond- 
ing exponent turns out to vary continuously with the 
mean percolation probability. The power-law behavior is 
due to 'blockages' at certain sites where the local perco- 
lation probability is small. Similarly, in the subcritical 
edge of the glassy phase, the spreading agent becomes 
localized at sites with high percolation probability. In 
d > 1, however, numerical simulations indicate that a 
glassy phase does not exist. 



B. Temporally quenched disorder 

Temporally quenched disorder is defined by the corre- 
lations 



Vt{t)vtit')=iS{t-t'). 



(All) 



In this case the additional term scales as a relevant per- 
turbation prjt -^ A^^/'^^^prjt- Therefore, we expect the 
critical behavior and the associated critical exponents to 
change entirely. In the field-theoretic formulation this 
corresponds to adding a term of the form 



S — Sn 



-7 



dt 



d xip^ 



(A12) 



The influence of spatio-temporally quenched disorder was 
investigated in detail by I. Jensen |7g]. Employing se- 
ries expansion techniques he demonstrated that the three 
exponents P^v^^^vn vary continuously with the disorder 
strength. Thus the transition no longer belongs to the 
DP universality class. A field-theoretic explanation of 
these findings is still missing. 



C. Spatio-temporally quenched disorder 

For spatio-temporally quenched disorder, the noise 
field 77 is uncorrelated in both space and time: 
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77.t(x, t)ry,t(x', t') = 7 d\x - x')<5(i - t') . 



(A13) 



In Rcggcon field theory, this would correspond to the 
addition of the term 



S = So+j d'^xdt 



(A14) 



which is an irrelevant perturbation. Spatio-temporally 
quenched disorder is expected in systems where each time 
step is associated with a new set of spatial degrees of free- 
dom. Examples include water in porous media subjected 
to a gravitational field as well as flowing sand on an in- 
clined plane. In these systems the critical behavior of DP 
should remain valid on large scales. 
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